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The relaxation function theory of doped two-dimensional 
5 = 1/2 Heisenberg antiferromagnetic (AF) system in the 
paramagnetic state is presented taking into account the hole 
subsystem as well as both the electron and AF correlations. 
The expression for fourth frequency moment of relaxation 
shape function is derived within the t — J model. The pre- 
sentation obeys rotational symmetry of the spin correlation 
functions and is valid for all wave vectors through the Bril- 
louin zone. The spin diffusion contribution to relaxation rates 
is evaluated and is shown to play a significant role in carrier 
free and doped antiferromagnet in agreement with exact diag- 
onalization calculations. At low temperatures the main con- 
tribution to the nuclear spin-lattice relaxation rate, ^^(l/Ti), 
of plane ®'^Cu arises from the AF fluctuations, and ^^(l/Ti), 
of plane ^^O, has the contributions from the wavevectors in 
the vicinity of (vr, tt) and small g ~ 0. It is shown that the the- 
ory is able to explain the main features of experimental data 
on temperature and doping dependence of in the 

paramagnetic state of both carrier free La2Cu04 and doped 
La2-a;Sr2,Cu04 compounds. 

PACS numbers: 74.72.Dn, 74.25.Ha, 75.40.Gb, 71.27.-ha 



I. INTRODUCTION 

The spin dynamics in doped two-dimensional S = 1/2 
Heisenberg antiferromagnetic (2DIIAF) systems remains 
the one of the intriguing problems of condensed matter in 
connection with physics of layered copper High Tempera- 
ture Superconductors (HTSC).^ The effect of doped holes 
in two-dimensional (2D) antiferromagnetic (AF) back- 
ground was studied in many papers.^ The Hohenberg- 
Mermin- Wagner theorem states the absence of long range 
order in low-dimensional isotropic Heisenberg systems 
at any finite temperature due to fluctuations, making 
the order short-ranged. The temperature dependence of 
correlation length in two-dimensional Heisenberg model 
was first described by Chakravarty, Halperin and Nel- 
son (CHN) by a quantum nonlinear a modeP in accord 
with neutron scattering (NS) experiments^ in carrier free 
La2Cu04. A significant advance was achieved in under- 
standing the 2D Heisenberg systems at low temperatures 
due to the improved further results of CHN in the renor- 
malized classical regime. Since then the temperature 
dependence of correlation length has been studied, e.g., 
by isotropic wave theory'' and by quantum consideration 
of skyrmions.* 

Nuclear Quadrupole Resonance (NQR) and Nuclear 
Magnetic Resonance (NMR) methods are very power- 
ful in studying the low energy excitations of HTSC.^ 
Chakravarty and Orbach^ developed a theory of magnetic 
relaxation phenomena in 2DHAF based on the quantum 



nonlinear a model of CHN in the critical region of fiuctu- 
ations. However, they considered only the contribution 
that arises from the wave vectors in the vicinity of AF 
wave vector (tt, tt). Despite the fact that contribution 
from the wave vectors q ~ does not dominate in the 
plane copper spin-lattice relaxation rate at low temper- 
atures and was not accounted in the analysis, the con- 
tribution from spin diffusion plays an important role in 
HTSC, especially in the plane oxygen relaxation rate as 
measured by NMR, since the AF fluctuating contribution 
is filtered out by the oxygen formfactor. 

The Nearly Antiferromagnetic Fermi Liquid (NAFL) 
model of Millis, Monien and Pines^^ reconciled the puz- 
zling observation of non-Korringa temperature depen- 
dence of the copper nuclear spin-lattice relaxation rate 
and Korringa temperature dependence of the oxygen 
and yttrium nuclear spin-lattice relaxation rates in op- 
timally doped YBa2Cu307 by postulating both the lo- 
calized Cu^+ magnetic moments and free oxygen holes. 
The NAFL model gave the relation between the AF cor- 
relation length and the relaxation rates and was applied 
in a wide temperature and doping range. However, this 
theory has some disadvantages, connected mainly with 
the phenomenological character of the NAFL description, 
since the temperature and doping dependence of corre- 
lation length was postulated or, at best, taken from a 
comparison with experiment. 

As a consequence, it is tempting to consider the 
2DHAF systems doped by charge carriers microscopically 
and the difficulties rise on. The numerically exact meth- 
ods study the relatively small clusters and in addition, 
these methods are hardly applicable if we need to obtain 
the dynamic quantities. In the absence of exact solution 
it is necessary to find a reliable approach that will de- 
scribe the physical quantities with convincing accuracy. 
This is intriguing especially since the observation of dras- 
tic change of various physical quantities with doping and 
emergence of "stripe" physics in doped HTSC and related 
compounds. 

The t — J model became very popular since it was 
pointed by Anderson^^ as a perspective model to describe 
the electronic properties of layered HTSC cuprates. The 
present theory is developed for paramagnetic state using 
a Mori-Zwanzig projection operator procedure, ^^'^^ with 
a three-pole approximation for the relaxation function. 
The merits of the theory were demonstrated^^ by com- 
parison with experiments in antiferromagnets in a wide 
temperature range down to temperatures close to Neel 
temperature T/y. The advantage of the present formula- 
tion is that it allows to take into account not only the 
AF correlation effects in relaxation rates, but also the 
contribution from spin diffusion. The presentation of the 
t — J model in terms of Hubbard operators is known 
to obey the rotational symmetry of the spin correlation 
functions and automatically guarantees the exclusion of 



double occupancy. The form of static susceptibility will 
be used from microscopic theory^^ as obtained beyond 
the Random Phase Approximation (RPA) and naturally 
takes into account the contribution from the hole subsys- 
tem. The dynamic structure factor is a quantity directly 
measured in various experiments, e.g., by magnetic neu- 
tron scattering. The advantages of expressing the dy- 
namic structure factor in terms of relaxation function 
were shown by Mori and Kawasaki.^** Until present the 
Mori-Zwanzig projection operator procedure was applied 
in connection with HTSC only to carrier free 2D 5 = 1/2 
AF system"'^^ with a three-pole approximation for the re- 
laxation function^^ and the correlation length was used 
from the results of CHN."^ It should be emphasized that 
the method developed in the present work for calcula- 
tions of dynamic structure factor and nuclear spin-lattice 
relaxation rates is similar, however, the approach is dif- 
ferent from the calculations with the approximations for 
dynamic spin susceptibilities. ^^'^'^^'^'''^^'^^ 

The paper is organized as follows. In Sec. II, the basic 
relations are presented for the relaxation function with 
a three pole approximation in a continued fraction rep- 
resentation of the Laplace transform and the static spin 
susceptibility. Sec. Ill shows the evaluation of the second 
and the fourth frequency moments of relaxation shape 
function within the t — J model. Sec. IV presents the re- 
sults of calculations, comparison with experiment, other 
theories and discussion. Sec. V is the Conclusion. 



II. BASIC RELATIONS 

We employ the t — J Hamiltonian written in terms of 
the Hubbard operators: 



1>J 



(1) 



Here, Si are spin-1/2 operators at the lattice sites i, and 
X^^ are the Hubbard operators that create an electron 
with spin a at site i. The hopping integral tjj describes 
the motion of electrons causing a change in their spins. 
In this paper, tij = t refers to hopping between nearest 
neighbors and J is the nearest-neighbor antifcrromag- 
netic (AF) coupling constant. The spin and density op- 
erators are defined as follows: 



Xf 



Hi 



[a 



(2) 



(3) 



with the standard normalization X^^ + X^"^-\-X~~ = 1. 
Without loss of generality, we measure all energies from 
the " center of gravity" of the band. 



A. Mori-Zwanzig projection operator procedure and 
a three pole approximation for the dynamic 
relaxation function 

In the present formulation we will follow Mori.^** The 
time evolution of a dynamical variable S^^{t), say, is given 
by 



dr 



(4) 



In general, L is the Liouville operator, and in our, quantal 
case, iLS^{T) is the corresponding commutator with the 

Hamiltonian (1). The projection of the vector S^jy) onto 
the = S^{t = 0) axis is given by 

PoS^{r) = R{k,T) ■ S^, (5) 

and defines the linear projection Hcrmitian operator Vq. 
One may separate S^{t) into the projective and vertical 
components with respect to the axis: 



S^{T) = RiKT)-S^ + {l-Po)S^iT), 



where 



(6) 



(7) 



is the relaxation function in the inner-product bracket 
notation: 

{S^{r),{S^_ur )^kBT dg 

X (eMQHW{T)eM-QH){Sl^r) , (8) 

where the angular brackets denote the thermal average. 

For future evaluations, it is convenient to introduce a 
set of quantities /o(t), /i(t), ... , fjir), . . . defined by 
equations 

/j(t) = ex.p{iLjT)fj = eKY>{iLjT)iLjfj_i, {j > 1), (9) 



where /o(t) 
and 



A? = •(/,•-!, /A)"'- 



(10) 



The set {fj} forms an orthogonal set. The larger number 
of fj is used, the finer description of S^{t) is obtained. 
The last quantity from this set /„, affected by evolution 
operator exp(zL„T), resulting in /^(t), was called the "n- 
th order random force", acting on the variable S^{t) 
and is responsible for fluctuation from its average motion. 

In terms of Laplace transform of the relaxation func- 
tion, i?^(k, t), one may construct a continued fraction 
representation for i?(k, s): 

POO 

i?^(k,s)= / rfTe^^i?(k,T) 

= l}{s + AlJ[s + AlJis + AlJ . . .)]}, (11) 
where A|]^ are related to the frequency moments 



J-oo * 



d"i?(k, r) 



(12) 



T=0 



of the relaxation shape function 



F{k,u!) = -Re[R^{k,iuj)] 



1 
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(13) 
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Lovesey and Meserve^^ truncated the relaxation func- 
tion (11) to third order. They argued that since the k 
dependence of A2k is much weaker than that of Aik, and 
using the analytical rcsults^"^ for the sixth frequency mo- 
ment (w^), the approximation of Aak by a constant is a 
good approximation. Thus, they suggested a three pole 
approximation for relaxation function, 

i?(k, s) = l/{s + AlJ[s + AlJis + 1/Tk)]}, (15) 



with a cutoff characteristic time 



Tk 



1/2 



For F{k,u!) this is equivalent to 



F(k,w) 



A2 
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(16) 



(17) 



Here one should note that F(k, lu) is real, normalized to 
unity dLjF{k, w) = 1 and even in both k and lj. 

The dynamic structure factor 5(k, oj) is related to 
the relaxation shape function F(k, w) through the 
fluctuation-dissipation theorem 



^(k,^) = 



27rwx(k) 



1 - exp{-(j/kBT) 



F{k,u;). 



(18) 



The only undefined quantity in the present formula- 
tion is the static spin susceptibility x(k) in (18). Until 
present this method was used^® to describe the paramag- 
netic state properties of (anti)ferromagnets with the form 
of static susceptibility, which was justified only at high 
temperatures (see also Refs. 24 and 19). In the present 
work we will employ the microscopic formula for static 
spin susceptibility^^ that is shown to work in the overall 
temperature range and properly takes into account the 
hole subsystem. 



B. Static spin susceptibility 

From Ref. 17, the expression for static spin suscepti- 
bility x(k) is straightforward, 



X(k) 



4|ci| 



J9-{9+ + 'yk) 



(19) 



and its structure is the same as in the isotropic spin- wave 
theory.^ The meaning of g+ is clear: it is related to the 
correlation length ^ via the expression 



(20) 



where a is a lattice unit. For carrier free AF system the 
relation (20) was obtained from the exponential decay 
of the spin-spin correlation function at large separations, 
whereas at finite doping the same expression was derived 
from the expansion of x(k) taken around the AF wave 
vector Q, however, taking now into account the contri- 
bution from the hole subsystem.^'' 



Cl 



7E(^^^^p)' C2 = ;2^E(^^'^^W)' (21) 



are the nearest and next-nearest neighbor spin correla- 
tion functions, respectively, the index p runs over nearest 
neighbors, and 

7k = - exp(ik/9) — i(cos k^a + cos kyu), (22) 



Q- = 4:az\ci\, z = 4 is the number of nearest neighbors 
for square lattice. The parameters a and (3 were intro- 
duced in the decoupling procedures for the higher-order 
Green's fimctions.^'"' The parameter a preserves the im- 
portant property that spin operators obey the relation 
(Sf)=3/4 which should hold at all temperatures. The 
numerical values for the temperature dependence of ^ 
were determined in carrier free La2Cu04 using J = 0.12 
eV and treating /3 as the only adjustable parameter. 
The best fit to experimental data, which were deduced 
from NS,* was obtained with f3 = 2.5. This value will be 
kept fixed in the present calculations. The original self- 
consistent theory of Kondo and Yamaji (KY)^^ with a 
— 13 = 1.705 fails in explanation of the absolute values of 
^. Since f3 enters in the combination f3c2, the increase of 
the values of next-nearest correlations causes the exten- 
sion of short-range AF order and hence the enhancement 
of ^ together with the spin stiffness constant ps-^"^ In 
T ^ limit, for both the carrier free and doped case, C2 
and g- are related as, 



1(1. 



12C2/3). 



(23) 



The reliability of the theory has been demonstrated by 
comparing the numerical values for ci, C2 and xs = 
X(k = 0) with Monte Carlo, Exact Diagonalization cal- 
culations and other theories. 

In the present calculations for small doping 6 we will 
use the expression for doping and temperature depen- 
dence of ^, given by,^'' 



kBT 



■eM'^nps/kBT). 



(24) 



As discussed in Ref. 17, ^ diverges at T = even in 
doped samples and this disagrees with experiment. This 
disagreement appears, probably, due to the overestima- 
tion of the role of AF correlations at low temperatures 
in the KY decoupling procedure. To mimic the low tem- 
perature behavior of the correlation length we will use 
the expression, as in Refs. 21, 22, resulting in effective 
correlation length ^eff , given by, 



(25) 



and 



Thus, the theory is able to explain the temperature and 
doping dependence of correlation length. The expression 
(25) is different from Keimer et al.^^ empirical equation, 
where ^ is given by the Hascnfratz-Nicdcrmaycr formula^ 
and hence, there is no influence of the hole subsystem on 
^. In contrast, in the present theory, ^ is affected by 
doped holes. Thus from now on we replace ^ by ^eff- 



C. Excitation spectrum 

The band evolution with doping remains the contro- 
versial topic. The usual parameter set in the t — J model 
is t = J/O.'d.^ The electronic and AF spin-spin corre- 
lation fmictions reduce the hoppings^^ resulting in ef- 
fective values. In the early proposal oi t — J model by 
Anderson^^'^® for description of properties of layered cop- 
per HTSC compounds the phenomenological relation was 
settled for the band width ~ 6t. Following the idea of 
Zhang and Rice^^ about copper-oxygen singlets forma- 
tion it was shown by Eremin et a/.,^" that it is possible 
to describe correctly the elementary excitations spectrum 
in cuprates. This singlet correlated band is analogous 
to upper Hubbard band with essential distinction - the 
subband splitting is much smaller compared to Hubbard 
model. Therefore it is possible to apply Hubbard for- 
malism without strict restriction on t and J values in 
t — J model. Taking into account the AF spin-spin cor- 
relations selfconsistently, resulting in effective hoppings, 
it was shown that the band width varies linearly with 
doping in a wide range from lightly doped to optimally 
doped compounds. ^"'^ To avoid confusion from approxi- 
mations the simple expression for effective hopping 



teff = 6 J/0.3, 



(26) 



will be employed in the following evaluations to match 
the insulator - metal transition. At high hole concentra- 
tions {S ~ 1), where the correlation effects are negligible, 
Eq. (26) gives, as it should, the value of band width for 
the noninter acting case. 
Thus, Ei^ is given by 



E]f, = 2teff{cos k^a + cos kya). 



(27) 



This expression resembles well also the values of singlet- 
correlated bandwidth « 0.4 eV (see Ref. 30) in optimally 
doped {6 « 0.15 per Cu site) layered copper HTSC as 
obtained from Angle Resolved Photoemission Electron 
Spectroscopy (ARPES).^^ 



III. EVALUATION OF EXPRESSIONS FOR 
FREQUENCY MOMENTS 



)/x(k). 



(29) 



The main effort necessary here is to obtain the expression 
for {i^^}- In this Section wo will start with evaluation of 
commutators and calculation of the thermodynamic av- 
erages, then introduce the decoupling procedures and, 
finally, present the result for {i^^}- The procedure will 
be tested by comparison of our expression for spin part 
with the existing result. The Xf'^ and operators 
are fermions and obey the anticommutation relations, 
whereas the S"' and X^'^ are bosonic-like and obey the 
commutation relations. The terms with the operators of 
different type are assumed to satisfy the commutation 
relations. The commutators with the products of oper- 
ators were decomposed on terms that contain commuta- 
tors and(or) anticommutators depending on the type of 
operators. 



A. Evaluation of commutators 

In order to calculate the second (w^) and the fourth 
(w^) frequency moments we first evaluate the commuta- 
tors 



and 



(30) 



(31) 



The commutator of expression given by Eq. (31) with 
the hopping term Ht in (1) is 



[['S'm. Ht] ^Ht] = ^ timO- 1 2 



■yaO \raa vOcr 



I vaO -yaa vOa 



xf (X™ 

ii„XfO(X^o + Xr)xr}, (32) 



whereas the commutator of expression given by Eq. (30) 
with the spin part Hj of (1) is 

w f ca Qz na ca na nz i nz aa ca aa qz na \ 

I ^ \ ^ T T fo C-^ Q-^ Q-^ 

"T g 2^ 'Jim'Jlm\'^>->i i->m>-'l ~ '^i '-'m'-'l ~ '^i '^m'-'l 

+ QZ Qa Qa Q(T QZ Q(7 i <~f(T Q(7 QZ QCT QZ Q(T \ / *J *J \ 

'^m'-'l "^i ~ '-'m'-'l '-'i + '-'i '-'m ~ '-'l '-'i '-'mJ- W"^/ 

The rest commutators of this type are as follows 



We now describe the procedure used to calculate the 
second (w^) ^^'^ fourth (w^) frequency moments of 
F{k,ij), by calculating directly the corresponding com- 
mutators in 



)/x(k), 



(28) 



+ {Jim - Jim) til {S^X^^Xi"' - Xi^X^^'S^) , (34) 



and 



(38) 



vaO vOtT vaa i vaO na vOa \raO vaa vOa 



-xrxf°x°r) 



"-(jO \rOcr \raa 



\^(tO vaa v'Oo' 



I \r \r[)(j C!<T \^U(7 I CCV^Q-V^Oct 

Xaa v'fO v'Oo' I vctO ca v'Oo' -yaO v'5'cr vOctn 



(35) 



We conclude this subsection with the emphasis, that our 
Eqs. (31)- (35) are still exact. Wc will restrict further 
ourselves and take into account the correlations between 
the first and the second neighbors only. The resulting 
form of calculated commutators suggests the types of 
necessary thermodynamic averages. The calculation of 
these thermodynamic averages together with the decou- 
pling procedures we need to employ in order to estimate 
the values of higher spin, transfer amplitude and density 
correlation functions will be presented in the following 
subsections. 



B. Thermodynamic averages 

To calculate the thermodynamic averages, we use the 
retarded Green's functions formalism. The equation of 
motion for a retarded Green's function {{A\B))i^ takes 
the form 



uj{{A\B))^ = {[A,BU) + {{[A,H]\B)l 



(36) 



where (...) denotes the thermal average. The standard re- 
lationship between correlation and Green's function may 
be written as 



(BA) = dcuf{u;){{A\B))^, 



(37) 



where f{uj) ~ [exp {uj /ksT) -\- 1]^^ is the Fermi function; 
the contour encircles the real axis without enclosing any 
poles of /(w). 

In general, Eq. (36) cannot be solved exactly and 
one needs some sort of approximation. To evaluate 
the Green's fimction {{[A,H]\B))^ in Eq. (36), one uses 
a decoupling scheme originally proposed by Roth"^^ for 
calculations on the Hubbard model. It can be shown 
that Roth's method is essentially equivalent to the Mori- 
Zwanzig projection technique"'"''"'^ and is strongly related 
to the moments method as applied to the evaluation of 
the spectral density of the Green's functions. "^^'"^^ Roth's 
method has been studied by many authors, ^'^'"^^'"^^ and 
became a general method to treat approximately the 
quasiparticle spectrum of an interacting system. The re- 
liability of the method has been demonstrated by com- 
parison with the exact diagonalization results. "^^ 

Roth's method"^^ implies that we seek a set of opera- 
tors An, which are believed to be the most relevant to 
describe the one-particle excitations of the system of in- 
terest. Also, it is assumed that, in some approximation, 
these operators obey the relations^^ 



where the parameters Knm are derived through a set of 
linear equations 



{[[An,H],Ai] + ) = Y,Knm{[Am,Af] + ). 



(39) 



Thus, it remains to define the operators An. Because, in 

the framework of the t — J model, the quasiparticles are 
described by the Hubbard operators , a set of oper- 
ators An contains only one operator A = . Hence, 
the matrix Knm is diagonal and also contains one ele- 
ment K = E^, where is the energy of of an electron 
with wave vector k and spin projection a. Consequently, 
Eqs. (38) and (39) become 



([[x°^i^],xn+) = ^;^([xo^x^°]+) 



(40) 



(41) 



In the 2D t — J model, long-range order is absent at any 
finite temperature and hence, E^ does not depend on a. 
Thus, we can replace E^ and E^ by i?k- 

For our evaluations we need the thermal averages of the 
following types: {X^^X^'') and {X^'^Xf"'). The aver- 
ages (spin-spin correlation functions) of the type (SfSi) 
were defined in the SubSec. II B and the calculation pro- 
cedure together with the numerical values will be outlined 
in the Sec. IV. 

First, one should note, that in the absence of long- 
range order, {X['^) does not depend on the site index 
and hence, according to Eq. (3), Tq = (Xf*) = {X^") = 
(1 - 6)/2 and Cq = {S^S^) = 1/4. 

The transfer amplitude between the first neighbors 
Ti = pli is given by 



crO yOct 



and may be calculated using the spectral theorem 

2^7k/k- 



-^1 5Z gk->^ 

k e~^^ +1 k 



(42) 



(43) 



The parameter 7i in Eq. (43) has been estimated in Ref. 
21, 



7, « 1 f 1 _ e--^) - 2^, ~5= ^ 



1 + 5' 



(44) 



with an accuracy of a few percent over the whole region 
of 6 from to 1 . Here one should note that for very small 
S and low temperatures, Ji w 2S. 

The transfer amplitude between the second neighbors. 



rp _ 1 \ ^ / vOct \ 



(45) 



is 



z{z-l) 



167^ — 4 cos k^a cos kytt — 4 



= - _ ^ (167^ - 4 cos fc^a cos kya-A)j^. 



For p we have 



1 + 5 

2 ' 



(46) 



(47) 



where 6 is the number of extra holes, due to doping, per 
one plane Cu^+. The chemical potential fi is related to 
<5by 

P_ 
N 



(48) 



where = [exp{—E]^ + ii)/kBT+l]~^ is the Fermi func- 
tion of holes. 

To obtain the thermodynamic averages of the type 
{Xf'^Xf_^p ) it is convenient to make the following def- 
initions 



and 



p 



(49) 



(50) 



To obtain A and X^s we use the two Green's functions^^ 

= ^E((^kn^f ^r.))- (51) 

p 

gL'H'^) = \Y.i{Xl'\^t'^i+p))-- (52) 
p 

Note that in the paramagnetic state, X^a = X^a and 

According to Eqs. (36) and (40), the equation of mo- 
tion for g'^^ (lu) and G^^ (to) can be written as 

{uj - Ek)G'^^ (a;) = - p - X^, + pljk) (53) 



\/]V 



(a; - Eu)g'^\co) = ^(1 - p - A^^). (54) 



\/]V 



According to Eq. (37): 



da;/(w)GW(w), (55) 



= ^ E ^ ^ da;/(a;)G;,'='(a;). (56) 

k 



Consequently, Eqs. (53) and (54) lead to a system of 
linear equations for A^g and A^cr with the trivial solution 

(57) 



x=x,, = {i-pr-^i' 



2p-l ' 



A<.g = (l-p-A)i-^ = (l-p)^ + 



2 , (1-P)P%2_ 

2p - 1 ^ ^ 



C. Decoupling procedures 

We now describe the decoupling procedures for the 
thermodynamic averages performed in spirit of Hubbard 
and Jain^^ and Kondo and Yamaji.^^ 

The averages of the type {XfX^^X^Xf) are decou- 
pled resulting in products of transfer amplitudes and the 
decoupling parameter C,, 

{XfX^''X~^Xf) ^ axfX^'^){X^^Xf). (59) 

The four-spin correlation functions are approxi- 
mated, as usually, by products of two-spin correlation 
functions, ""^^ however, multiplied now with the decoupling 
parameter C,. Thus, we employ the decoupling procedures 



and 



(60) 



(61) 



for i r and m ^ j, whereas 

{S;:S^S:^SJ) ^ 2co{S^^Sf). (62) 

The averages with the products of operators Xf^X^"^ 
between the nearest (next-nearest) neighbors with (1 — 
X^^){1 - Xf"') are decoupled as follows: 

(XfOXO-(l-X^-)(l-X/-')) 

^ (XfOXO-)(l - - Xf-' + X^JXf-'). (63) 
and so on. 

The averages with spin and Hubbard operators are de- 
coupled as follows: 

{XfX'^j^SfS"^) ^ (Xf"xO'^)(5f5;), (64) 

and with spin and density operators: 



I ~v-cc QU Qrj\ ^ / f \ / CO" Q<J\ 



(65) 



The averages {X^'^X'^ ) between the second neigh- 
boring operators are decoupled simply by 



(xrxf+5)-(xr)(xf4), 



(66) 



because an inspection of Eqs. (57) and (58) shows that 

the values of averages of these type between the first 
neighbors differ only slightly from {X"") {X"^^^. There- 
fore, the averages between second neighbors in Eq. (66) 
are thought as independent. In addition, because the 
averages (Xf'^Xj' ) between the first, in contrast to 
next-nearest, neighbors, are calculated exactly, the av- 
erages like {X^^X^X" " ) are decoupled in a way to 
avoid, where possible, the averages of the type as given 
in Eq. (66). 



D. Final result for (wk) and (oj^) 

The expression for the second moment (w^) is straight- 
forward. Calculating the commutator with the expres- 
sions given by (30) and (31) with S^, taking the thermal 
average and the Fourier transform, the result is 

(col) = - (8Jci - itejfph) (1 - 7k) /x(k). (67) 

We now proceed with calculating the commutators 
for (t^k)- Taking the commutators with the expressions 
given by Eqs. (32)-(35) with Eqs. (30) and (31) wc ob- 
tain the expression for (w^). The expression was then 
approximated using the decoupling procedures for ther- 
mal averages as was described in the previous SubScction. 
Taking the Fourier transform we arrive to expression that 
contains various types of sums over lattice sites. The cor- 
responding types of sums over lattice sites are presented 
in Appendix A. Utilizing these sums we obtain 

(o;^) ^ -{128J3[C2 (1 - 7^) (CC2 (7k - |) - \co) 
+ coci (I - i7k + §7^) + CciC2 - f 7k + ^7^) 
+ Cc? (1 - f 7k + x7k + I cos k^a cos kyU - 2-fl)] 
+ I6plitlff[ci (3 — 27^ — coskxacoskyO) 
+ CT2 (7 - 127k + 57^) + ^ (1 - 47k + 37^) 
+ {S + X) (-^ + 97k - 37k ~ 1 cos k^a cos kyo) ] 
+ m^ffJ'ph[co (-f + f7k-f7^) 
-|- ci (l67k ~ "^^^k + 257k ~ I ~ I cos k^a cos kyu) 
+ C2(-f + f7k-^7^) + ^i^(7k-l)] 

+ mtlffJic (|7k - !) + (67^ - f 7k + f ) 

+ T2 (iTo + |r2 - A) (272 - 37k + 1) 



+ 



literature. One should note that [S^,ni] = 0. Thus, 
the commutator (30) does not depend on the density- 
density term riirij in (1). Moreover, the commutator of 
spin operators appearing in (30) with the density 
operators rn is zero also: [S!^,ni] = 0. This enables us 
to compare the result (68) for spin part with the results 
for carrier free 2DHAF. As it was already mentioned, we 
restrict ourselves and take into account the correlations 
between the first and the next-nearest neighbors only. 

For the second moment {ui^) the result (67) is evi- 
dently correct. To compare the results for (t^^) we need 
some additional effort. The corresponding expressions 
and the procedure will be described in Appendix B. The 
expression, as shown in (68), with teff = and with the 
decoupling parameter C settled to unity C = 1 is the same, 
as given in Ref. 16 and may be compared using the inte- 
grals (sums) over the Brillouin zone (see Appendix B). 

IV. COMPARISON WITH EXPERIMENT AND 
DISCUSSION 

The results of the calculations are summarized in the 
Table 1. The value of extra holes, due to doping, per one 
plane Cu^^, S, can be identified with the Sr content x in 
La2-a;Srj;Cu04. The AF spin-spin correlation fimctions 
Ci, C2, the spin stiffness constant ps and the parameter g- 
were calculated using the expressions and the procedure 
as described in Ref. 17 (see also Rcf. 21) in the T ^ 
limit, since wc will employ them in a temperature range 
T < 1000 K ~ 0.7 J, where experimental data exist and 
according to the calculations their values have a weak 
temperature dependence. 

Table 1. The calculated at T ^ values of AF spin- 
spin correlation function ci between the first neighbors, 
the g_ parameter and the spin stiffness constant ps as 
a function of La2-a;Sra;Cu04 doping x together with the 
values of decoupling parameter ( as extracted from com- 
parison with ^^Cu NQR spin-lattice relaxation rate mea- 
surements. 



X ci 2-!Tps/ J C 



+ (f A+- (1 - To) - T2 + ATo) (7k - 1) 
+ (27^ - |7k + I) + c,To (f 7^ - i7k + i) 

+ CiT2(^7^-f7k + f) 
+ C2To(-272 + |7k-|) 
C(p/i)^ (-47k + 67k + T7k - cosk^acoskya - ^) 

+ T2C2(l67^-2l7^-|7k + ¥) 
+ T0T2 (27^ - |7k + I) + T2C0 (-57^ + |7k + i) 

+ CT| (-27^ + 67^ - f 7k + |)]}/X(k). (68) 
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To insure the accuracy of the calculations we compare 
the result for spin part Hj with that reported in the 



A. Temperature and doping dependence of 
antiferromagnetic correlation length 

The AF correlation length, its doping and tempera- 
ture dependence is given by Eq. (25). is the value 
of correlation length at T = and gives information on 
the topology of holes. At present, two types of doping 
dependencies are under debate in the literature. The 
first, *° is the localization of holes close to the randomly 
distributed Sr^+ ions gives = aj \fx.. The other one 
is the formation of dynamic domain walls. In this case 
^0 = a/nx, where n is the average distance between the 
holes along the domain walls in lattice units,**^'**^ usually 
called "stripes". 

The best fit of i^^ff to experimental data^^ yields = 
a/nx, where n = 1.3 for samples with x ^ 0.02 and n = 2 
for samples with x > 0.02, see Figure 1. These results 
arc in agreement with the the conclusion of Borsa et al.'^^ 
for La2-2;Sr2;Cu04 compounds with x ^ 0.02 and with 
the analysis of Carretta et al. on the basis of dynamical 
scaling.**'^ The presence of stripes with n = 2 was found 
also by Tranquada et al. in a; ~ 1/8 compound. 
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FIG. 1. The inverse correlation length versus tem- 
perature fitted (solid lines) to the experimental data as ob- 
tained from neutron scattering experiments. For carrier free 
La2Cu04: filled circles from Ref. 4 (the fitted data), aster- 
isks from Ref. 44 and open circles from Ref. 26. For doped 
La2-xSra,Cu04: squares for x = 0.02, down triangles for 
X = 0.03 and up triangles for x = 0.04 from Ref. 26. 

The fit with the relation a/nx in may be obtained 
by comparing the experimental data in the low tem- 
perature region. Indeed, the doping relation a/ \/x fits 
with the neutron scattering data poorly than the a/nx 
does, keeping in mind the value of plane lattice con- 
stant a = 3.79 A. For example, for x = 0.02 one has 
^0^^ = v^/3.79 = 0.037 A,~^ whereas the experimental 
value is 0.005-0.007 A.-^ The inability to fit the l/^/x 
relation within the NS data^^ is clear also for larger val- 
ues of doping: ^^"^ = V^/3.79 = 0.046 A^^ for x = 0.03 
and ^(7^ = V^/3.79 = 0.053 A^^ for x = 0.04 differs 
strongly from 0.013-0.02 A"^ and 0.02-0.026 A,-\ con- 
sequently, as obtained by NS^^ (see Fig. 1). On the other 
hand, the relation — nx/a gives good agreement: 
0.0069 A-i for X = 0.02 with n = 1.3, and 0.016 A"! 



for X = 0.03, and 0.021 A^^ for x = 0.04 with n = 2. 
The value n = 1.3 for x = 0.02 fits the values of corre- 
lation length as obtained from experimental data better 
in contrast to that with n = 2. This result is especially 
evident when one tries to compare the values of at 
high temperatures. The high quality of the fit is in agree- 
ment with previous studies and seems to confirm the mi- 
croscgrcgation, however, with different values of average 
distance between the holes along the domain walls. It 
is tempting to speculate that the change in the average 
distance n from n w 1.3 for a; < Xc and n = 2 for a; > Xc 
appears at the value of Sr content when the Nccl order 
is completely suppressed (T^r ^ at x = Xc — 0.02) as 
discussed, consequently, for example, in Refs. 41 and 45. 

Having established the temperature and doping depen- 
dence of correlation length we now proceed with calcula- 
tions of the nuclear spin-lattice relaxation rate. 



B. Spin-diffusion constant 

The spin-diffusion is described by the (k,aj) = (0,0) 
mode. The spin- diffusion plays an important role in an 
antiferromagnet and makes a pronounced contribution 
to the relaxation rates. The spin-diflFusion constant D is 
given by^^ 



D = lim 



limlJ?M)!. (69) 



k^o 7rfc2F(k, 0) k^o A;2 y 2 (w^) ' 
For small q, Af^^g tends to zero as 

A?q_.o = - + 9'«')(2Jci - teffpli)/xs, (70) 
whereas A^^^o — ('^q^o) / ("^q^o) remains finite, 



Aiq^o = { 



128J^ 



hc2{CC2 - Co) - jCciC2 + jCoCi 



1 1-5 

2 2 



+ 16iejJVl(|c2-ico-^i-^) 



+ 16tlffJ 



-AiMci-iT2(iTo + fTo2-A) 



- i(|A+-(l - To) - T2 + ATo) + ic2i±^ - ^ci 

- iciTo - ic2To - ^({phf + iciTa - IT2C2 

+ \T2T0 + ^T2Co - ^CT|] } /{2Jci - t^ffph). (71) 

The values of D for carrier free AF system may be 
compared with the results of other theories. The infi- 
nite temperature result is the same as in Ref. 19, namely 
D{T 00) = lV2nJa'^ = O.SlSSJa^. This value 
is close to that obtained by Morita,*^ D = 0.43Ja^ 
and D = ^V^nJa"^ = 0.3963 Ja^ in Refs. 23, 48, 49. 
The calculated value of £> = 2A6Ja'^ at T ^ (with 
^ = 1) is larger compared to that obtained in Ref. 19, 
D = 1.63 Ja2. The value of I? = 2.66Ja2 (with C = 1-8) 
is compatible with D w 3 Ja^ obtained at T = 900 K in 



Rcf. 48. In general, D weakly changes with doping and 
when one vary Figure 2 shows the doping dependence 
of D for two cases, with Q = 1 and when Q is obtained 
from the best fit to NQR data. The calculated doping de- 
pendence of D in the T — * limit may be compared with 
the results of Bonca and Jaklic^° at high temperatures, 
assuming that the doping dependence of D remains the 
same. Indeed, Figure 2 shows the remarkable agreement. 
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FIG. 2. The calculated doping dependence of the spin dif- 
fusion constant D in the T — » limit with the fixed decou- 
pling parameter C = 1 (squares) and with C, as obtained from 
the best fit to NQR data (circles). The solid line is the re- 
sult of Bonca and Jaklic^" for doping dependence of D in the 
high temperature limit. The asterisk marks the value of the 
spin diffusion constant D in the limit of high temperatures 
T ^ oo. 



C. Dynamic structure factor 

Now, since the nuclear spin-lattice relaxation rate is 

determined by the dynamic structure factor S{(\^uj) we 
first consider its temperature, frequency and wavevector 
dependence. First, we note that for all temperatures the 
relaxation shape function F(q, w) as well as S'(q, w) give 
the elastic peak at g = and u = Q. This is clear from 
Eq. (17) since both (ci^q) and (wq) ^ for small q. 
We now turn to the case of finite but small NMR/NQR 
frequencies. 

The relaxation shape function F{]<i,uj) with small lo, 
compared to temperature scale of the system, related to 
5(k, w) as 



S'(k,w ~ 0) = 27rfcBTx(k)F(k,c 



(72) 



We now explore the form of 5'(k, w) and obtain two peaks: 
one at q ~ and the other at Q = (7r/a,7r/a). The 
5(q, u) value at small q is large when the q values are 
such that 



A?q ^ a;TqA2q, 



(73) 



see Eqs. (17), (18), and (70), (71). Thus, ^(q.a;) has a 
sharp peak at qo: 



^(qo,^) 



ksTxs 



(74) 



and the 5(k, w) value at Q = {w/ajiT/a) is given by 



S{Q,Lo) 



2fcsTx(Q)TQA2Q 



^IQ 



(75) 



The relation given in Eq. (74) is in agreement with the re- 
sult of Makivic and Jarrel^^ on frequency dependence of 
the dynamic structure factor at small values of wavevec- 
tors as extracted from combination of the Maximum En- 
tropy Method and Quantum Monte Carlo calculations. 
This result agrees also with the basic relations known in 
the literature. From general physical grounds, namely, 
linear response theory, hydrodynamics and fluctuation- 
dissipation theorem, the diffusive spin dynamics gives the 
form of dynamic structure factor^^ 



5(q~ 0,w ~ 0) 



2X6 



1 - exp(-a;/fcsT) to"^ + {Dq^y 



(76) 



for small q and uj. 

Using Eq. (76) (or, equivalently, from Eq. (73)) one 
may easy estimate the value of qo which is given by 



4 



(77) 



For typical value of the measuring frequency, w « 1 mK, 
qoa ^ IT xlO~^ and weakly changes when one vary doping 
and the decoupling parameter C ((jon = qoy = "^^o)- 

For small q <C qo with finite ui the relaxation 
shape function F(q, w) and the dynamic structure fac- 
tor 5(q, a;) approaches zero: S'(qo ^ q ^ 0,cj) ^ 0. 

Thus, the contribution to the nuclear relaxation rates 
from q around has no peculiarities since w is finite, 
but small, compared to any values of the variables (we 
use J = 0.12 eV (1393 K)'). Figures (3)-(10) show the 
calculated dynamic structure factor for different values 
of doping, temperature and decoupling parameter ( with 
Wc = 27rx 33 MHz ( = 1.365 x 10"'^ eV) and with ujq = 
27rx 52 MHz ( = 2.15 x 10"^ eV) for x = 0.035. 
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FIG. 3. Log-scale mesh of the calculated dynamic structure 
factor ^(k, cjc) for carrier free antifcrromagnet at T = 500 K 
with simple decoupling C = 1- The cross on the vertical axis 
marks the maximum of JS'(k,Wc) at k = qo. 
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FIG. 4. Log-scale mesh of the calculated dynamic structure 
factor S(k,a;c) for carrier free aritifcrrornagnet at T = 1000 K 
with simple decoupling C = 1- The cross on the vertical axis 
marks the maximum of JS{k.,(jOc) at k = qo. 
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FIG. 5. Semilog-scale mesh of the calculated dynamic 
structure factor 5'(k,a;c) for carrier free antiferromagnet at 
T = 500 K with the decoupling parameter C, = 1.8. The region 
of small q values (q ^ qo) for which Logio[J5'(k, ui)] < —3 is 
not shown. The cross on the vertical axis marks the maximum 
of J5(k,a;c) at k = qo. 
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FIG. 6. Semilog-scale mesh of the calculated dynamic 
structure factor S{k.,ijJc) for carrier free antiferromagnet at 
T — 1000 K with the decoupling parameter — 1.8. 
The region of small q values (q -C qo) for which 
Logic [JS'(k,aj)] < —3 is not shown. The cross on the ver- 
tical axis marks the maximum of JyS'(k, Wc) at k = qo. 
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FIG. 7. Log-scale mesh of the calculated dynamic structure 
factor ^(k, cuc) for x = 0.04 hole content at T = 500 K with 
simple decoupling (" = 1- The cross on the vertical axis marks 
the maximum of JS(k., uic) at k = qo. 
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FIG. 8. Log-scale mesh of the calculated dynamic structure 
factor 5(k,Wc) for x = 0.04 hole content at T = 1000 K with 
simple decoupling C = 1- The cross on the vertical axis marks 
the maximum of J5'(k, Wc) at k = qo. 
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FIG. 9. Log-scale mesh of the calculated dynamic struc- 
ture factor S'(k,iJo) for x = 0.035 hole content at T/J = 0.1 
(T ~ 140 K) with simple decoupling C = 1- The cross on the 
vertical axis marks the maximum of JyS'(k, wo) at k = qo. 
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FIG. 10. Semilog-scale mesh of the calculated dynamic 
structure factor S(k., ujo) for x = 0.035 hole content at 
T/J = 0.2 (T ~ 280 K). The region of small q values (q < qo) 
for which LogioiJS(k.,uj)) < —3 is not shown. The cross on 
the vertical axis marks the maximum of JyS(k, uio) at k = qo. 



D. Nuclear spin-lattice relaxation 

The nuclear spin-lattice relaxation rate "(1/Ti) is 
given by 



'(l/Ti) = -$]"nk)'5(k,cu), 



(78) 



where oj is the measuring NMR/NQR frequency, "_F(k) 
is the wave vector dependent hyperfine formfactor^^'^^ 



17 



F{kf = 2C^ (l + 7k), 



for planar ^''O and 



63 



F(k)2 = (^„(, + 47kS)2, 



(79) 



(80) 



for ^"^Cu sites. Aat and B are the Cu on-site and trans- 
ferred hyperfine couplings, respectively. The quantiza- 
tion axis of the electric field gradient coincides with 



the crystal axis c which is perpendicular to Cu02 
planes defined by a and b. For C in Eq. (79), we 
adopted the formula = g(C| + C^), where Cy 
and Cx are the plane oxygen hyperfine couplings for 
two axis perpendicular to c, and that includes the 
factor I to account ^''(l/Ti) as measured by NMR. 
The values of hyperfine couplings were taken as fol- 
lows: Aab = 3.7x lO"'^ eV, for transferred hyperfine cou- 
pling B the relation B = (1-^2.75 x)x 3.8x10"'^ eV is 
used to match the weak changes with Sr doping^^ and 
C = 2.8x 10~^ eV in accord with the values as extracted 
from NMR data and used in calculations of relaxation 

rates l.H. 19.43,45, 48, 53-56 

The ^''O and ^^Cu nuclear relaxation rates are essen- 
tially determined by the corresponding formfactors given 
by Eqs. (79) and (80). Figures 11 and 12 show the quan- 
tity "i^(k)^S'(k, w) in a log-scale plot versus {ky — 0) 
and a semilog-scale plot along the diagonal of the Bril- 
louin zone. The calculations show the 1/k^ wavevector 
dependence of S'(k, oj) for k > qo and the k^ dependence 
of ^(k,^) for k < qo. The form of "F(k)2S'(k, w) gives 
the peaks at qo and k ~ (7r/a, tt/o). Thus, two types of 
contributions dominate in the nuclear spin-lattice relax- 
ation rates. 

Here one should note, that the theory contains the de- 
coupling procedures for correlation functions and the de- 
coupling parameter has been introduced. In general, 
despite the complicated structure of (i^^), C allows to reg- 
ulate the contributions with k ~ (7r/a,7r/a) to the spin- 
lattice relaxation rate 1/Ti. The larger (, the smaller is 
the contribution from k ~ (7r/a,7r/o) to l/Ti. 

The contribution to 1/Ti from q around has no 
peculiarities because the measuring frequency is finite 
(Wc « 27rx 33 MHz ( = 1.58 mK) in ^'^Cu NQR and 
uo w 27rx 52 MHz ( = 2.5 mK) in ^'^O NMR measure- 
ments), but small, compared to any values of the vari- 
ables (we take J=0.12 eV (1393 K)). The calculations 
also show that "F(k)^S'(k, uj) has a broad local minimum 
at cim which values are given by + ~ tt'^/o^. 

A direct numerical integration over k is difficult, be- 
cause "i^(k)^S'(k, uj) has an extremely sharp peak at very 
small qo- This requires an unattainably large number of 
points in numerical integration over Brillouin zone. 

We first estimate the value of contribution to 1/Ti 
from small q because Figures (11), (12) show that the 
pronounced contribution to relaxation rates should come 
from excitations with q 0. In the formulation of 
Ref. 10, the diffusive contribution. 



(81) 



q~0 



is divergent. While usually, this logarithmic divergence 
of {l/Ti)oiff was argued to cut ofF,^"''** since in a real 
systems the mechanism destroying the diffusion should 
be taken into account. In La2_a;Sra;Cu04 the combina- 
tion of three-dimensional effects, finite length scale and 
the presence of disorder in Cu02 planes may suppress 
the effect of spin diffusion on spin-lattice relaxation rates. 
For example, using exact diagonalization technique in the 
high temperature region, it was argued that strong local 
perturbation induced by oxygen defects limit the spin- 
diffusion.^* Despite of considerable effort, the appropri- 



ate account of these effects is still unavailable because of 
lack of any exact analytical result. 

In the present theory the meaning of go is clear: it is 
the value of k at which the maximum in S{k, co) occurs 
and for small u) it may be treated as the cut off in the 
integration for calculation of the contribution from spin 
diffusion (see Eq. (81) and Ref. 10). For q <C qo the inte- 
gration over q should be taken over the factor ~ q^/q§. 
One should note, that in the present theory the form of 
the calculated dynamic structure factor S{k,ui) comes 
not from, e.g., disorder, but rather than from the relax- 
ation function and is in agreement with linear response 
theory, hydrodynamics and fluctuation-dissipation theo- 
rem and this result does not depend on the order of the 
pole approximation in the relaxation function theory. 

Expanding ^(q, lo) around go we obtain, 



{l/TijDiff ^ —\n{const x D/lj). 



(84) 



_ -F{0rkBTa^Xs , 



(82) 



where A depends on frequency through go. A simple and 
rough estimate gives 



A ~ ln(l/go) ln{const x J/ui). 



(83) 



This result explains the reason why the oxygen ^'^(1/Ti) 
relaxation rate as measured by NMR remains unchanged 
at 9 T {too = 277 X 52 MHz) and 14.1 T (wo = 
27rx 81.4 MHz) within the experimental accuracy.^^ One 
should note that oj is much less than J = 1.8x10*^ MHz, 
hence ln(J/52 MHz)/ ln( J/81.4 MHz) « 1.03. Indeed, 
a sophisticated calculation gives A (33 MHz) = 2.52, 
A(52 MHz) = 2.44, and A(81.4 MHz) = 2.37 and its 
value changes on less than 1% when one vary and dop- 
ing within X ^ 0.04. In view of the result that the spin 
diffusion contribution is 70 %, the relative shift of the 
measured ^^(1/Ti) will be « 2 %, that lies within the 
experimental error (see Figure 13). 

The ^''(1/Ti) relaxation rate has a contribution due an- 
tiferromagnetic correlations between copper spins, how- 
ever, for wave vectors only in the vicinity of Q = 
(tt/o, 7r/a) because the formfactor ^'^F{k) filters out the 
contributions with Q = (tt/a, 7r/a). This filter causes 
also the minor sensitivity to the decoupling parameter 
C- The value of ( = 1-3 for x = 0.035 (see Table 1) is 
a plausible guess. This contribution was calculated by 
direct summation over k in the region k > q^. 

Figure 13 shows the calculated ^^O relaxation rate in 
the sector of lightly damped spin waves at low temper- 
atures. It is seen that ^^O relaxation rate has a weak 
frequency dependence. The good agreement of calcu- 
lations (even without adjusting the parameter C) with 
experiment shows, that in case the mechanism destroy- 
ing the diffusion is present in La2-xSrj,Cu04, whether it 
is caused by the three-dimensional effects, finite length 
scale or the presence of disorder in Cu02 planes, it seems 
that it affects the spin-lattice relaxation rates only little 
through the contribution from spin diffusion. 

A fair agreement between experiment and calcu- 
lated ^^(l/Ti) for hghtly doped La2-,,:Sra;Cu04 may be 
thought as fortuitous because of the following reason. In 
the present theory, the contribution to relaxation rate 
from small q depends on spin diffusion constant D as, 



In the T — !■ limit the spin diffusion constant D 
has to diverge for both carrier free and lightly doped 
La2-a;Sra;Cu04 since both Heisenberg and t — J models 
have nonzero spin stiffness. ^'^ On the other hand, in the 
quantum critical region, ps < T < J (note that in the 
present theory ps — 0.06 J « 80 K and weakly decreases 
with light doping), the spin-diffusion constant scales as 
D - (see Ref. 3). For hghtly doped x w 0.03 sys- 

tems one may expect this scaling to be valid also, but 
with finite correlation length for T ^ ps (see Figure 1). 
This justifies the validity of the present formulation and 
the propriety of the results obtained for finite doping and 
low temperatures. 

It is worth to mention the results of calculations by 
Chakravarty et al.^'^^'^^ of relaxation rates for La2Cu04 
on various nuclei. The results of CHN gave the tempera- 
ture dependence of correlation length with the two-loop 
corrections in perfect agreement with experiment. The 
calculated plane copper relaxation rate ^"^{l/Ti) within 
the CHN theory is in agreement with experiment,'"'* how- 
ever, for plane oxygen the results of calculations does not 
reproduce the values of the measured ^^(1/Ti) for insulat- 
ing Sr2Cu02Cl2 (see Ref. 57). Moreover, the agreement 
between theory and experiment becomes worse when the 
two- loop correction has been taken into account. ^° A pos- 
sible reason of this inconsistency is the missed contribu- 
tion from spin diffusion^*^ and the electronic structure of 
Sr2Cu02Cl2 with pockets centered around (7r/2,7r/2) as 
observed by LaRosa et al.^^ Thus, for Sr2Cu02Cl2 one 
has to take into account the hoppings between the next 
nearest and next-next nearest neighbors that are beyond 
our present consideration. 

We now discuss the temperature and doping depen- 
dence of plane copper nuclear spin-lattice relaxation rate 
Figure 14 shows the fitted to experi- 

mental data in hghtly doped La2-a;Sra,Cu04. It is seen 
that the agreement is good for temperatures T < J/2. 
For temperatures T > J/2 the agreement is less satisfac- 
tory. This is a bit puzzling but could be due to better 
account of the ^ values in the low temperature region 
(T < J/2) compared to T ss J and the preexponential 
factor ^ 1/T in (24) and in the spin-wave theory,'' which 
is an artifact of the mean-field approach. For carrier free 
La2Cu04 and at T = 1000 K (T ~ 0.7J), C « 5a. Ob- 
viously, the value of correlation length ^, as extracted 
in the limit of large separations^^ satisfies the inequal- 
ity C^/o^ 1. At temperatures T < J the behavior of 
^^(1/Ti) is smooth: the slope of the curve decreases with 
increasing T. At high temperatures, ^■'(1/ri) is expected 
to have a minimum at T r:! J and to be dominated by 
spin diffusion at T > J. The validity of Eq. (24) does 
not allow to calculate l/Ti at T > J. Similar temper- 
ature dependence of ^^(1/Ti) was found in Ref. 19. It 
was found that the saturation of ^'^(1/Ti) at T ~ J/2 is 
"primarily due to a competition from the spin diffusion 
over the critical slowing down" . The present work shows 
that the overwhelming contribution arises from strong 
short-range antiferromagnetic correlations between cop- 
per spins and at low temperatures the contribution from 
spin diffusion to ®^(1/Ti) is small (see Figure 14). 



One may wonder about the results of calculations if 
we avoid fitting Fi ffure f 5 shows the results of 

calculations for with fixed decoupling parameter 

^ = 1. It is seen, that the theory is able to reproduce the 
main features of doping and temperature dependence of 
however, in poor agreement with respect to the 
numerical values. 

It is interesting to compare the results of the calcula- 
tions for carrier free AF system with other theories. After 
the seminal work,"^ Chakravarty and Orbach^ predicted 
the decrease of with increasing temperature at 

low temperatures as ^'^{l/Ti) - T^/^ cxp{2TTps/kBT). 
After passing through a wide minimum at T ^ J/2 
the calculated increases with T at high tem- 

peratures where the system is recognized to be in the 
quantum critical (QC) region. On the other hand, using 
the 1 /N expansion method on the N-component nonlin- 
ear sigma model an apparent formula was obtained for 
^^(l/Ti) in the QC region.^ However, the behavior of 
^^(l/Ti) was found to be nearly independent on tem- 
perature: ~ (T/J)0 °28. The calculations with 
thermally excited skyrmions by Belov and Kochelaev^^ 
showed that ^^(l/Ti) has also the nonmonotonic temper- 
ature dependence with the slight increase at T > 0.67J. 
In the present theory the calculated value of contribu- 
tion to ^^(1/Ti) relaxation rate from spin diffusion at 
T = 1000 K is 550 sec~^ in perfect agreement with 300- 
600 sec~-^ as estimated by Sokol et al^^ using exact diag- 
onalization technique. However, their estimate was based 
on the presence of disorder in Cu02 planes. 
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FIG. 11. Log-scale plot of the quantity °'F{i!.fS(k,u>)/J 
versus {ky = 0). The solid lines from down till up: 
^'^F(kfS(k,ujo)/J is given for x = 0.035 at T/J = 0.1 
(T ~ 140 K) and T/J = 0.2 (T ~ 280 K) with the de- 
coupling parameter ^ = 1.3, ®^F(k)^S(k, Wc)/J is given for 
carrier free AF (a; = 0) at T = 500 K and T = 1000 K 
with the decoupling parameter — 1.8. The dotted 
lines are "F(q())^S'(q(), a;)fci:^/(0.5Jqo) for kx < qo and 
"F(qo)"S(qo,w)qi/(0.5Jfe.") for > qo. 
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FIG. 12. Scmilog-scale 

plot of the quantity ™F(k)^S(k, cj)/ J versus k along the di- 
agonal of the Brillouin zone (fc, = ky). i^F(k)2S(k,wo)/J 
is given for x = 0.035 at T/J = 0.1 (T ~ 140 K, lower 
curve) and T/J = 0.2 (T ~ 280 K, upper curve) with the 
decoupling parameter ( — 1.3. ^^-F(k)^S(k, (j;c)/J is given 
for carrier free AF (a; = 0) at T = 500 K (lower curve at 
small k) and T — 1000 K (upper curve at small k) with 
the decoupling parameter ^ = 1.8. The dotted lines are 
"F(qo)^S'(qo, ti;)qo/(0.5 Jk^) for k > go. The lower and upper 
crosses mark the values of ''•^F(qo)^S'(qo, a;c)/J at T = 500 K 
and T = 1000 K, respectively. 
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FIG. 13. The calculated temperature and doping depen- 
dence of the plane oxygen nuclear spin-lattice relaxation rate 
^^{1/Ti) (lines) and the experimental data for La2-a;Sra^Cu04 
as measured by NMR with x = 0.025 (triangles) and 
X = 0.035 (squares) from Ref. 57. The experimental points 
have been rearranged with J = 1393 K. The results of calcula- 
tions with uj = 27rx 52 MHz (9 Tesla) are given for x = 0.035 
(C = 1.3) by solid line and for x = 0.025 (C = 1.6) by dashed 
line. The result of calculation with ui — 2nx 81.4 MHz for 
X = 0.035 = 1.3) coincides with the dashed line. The con- 
tribution to ^^(1/Ti) from the wave vectors in the vicinity of 
(7r/a,7r/o) for x = 0.035 with = 1 is shown by upper dotted 
line and by lower dotted line with ( = 1.3. 
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FIG. 14. Temperature and doping dependence of the plane 
copper nuclear spin-lattice relaxation rate ''"^(l/Ti) — 2W 
(solid lines) fitted to the experimental data for La2-:ESr:ECu04 
with x = 0, x = 0.012, X = 0.024 and x = 0.03 from Ref. 45 
and with x = 0.04 from Ref. 58. The values of the fitting 
parameter are shown in the Table 1. The dashed line shows 
the contribution to ''^(1/Ti) from spin diffusion for a: = 0. 
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FIG. 15. Recalculated plane copper nuclear spin-lattice re- 
laxation rate ^^(1/Ti) = 2W with fixed ^ = 1 is shown by 
solid lines from up till down with increasing doping. The 
dashed line shows the contribution to from spin dif- 

fusion for X = 0.04 with C = 1- The notation is the same as 
in Fig. 14. 



n ^ 1.3 for samples with x ^ 0.02 and n = 2 for samples 
with X > 0.02. The expression for fourth frequency mo- 
ment of relaxation shape function is derived within the 
t—J model. The spin diffusion contribution to relaxation 
rates is evaluated and is shown to play a significant role 
in carrier free and doped antiferromagnet in agreement 
with exact diagonalization calculations. The convergence 
of contribution from spin diffusion to spin-lattice relax- 
ation rates is preserved by linear response theory and 
hydrodynamics. At low temperatures the nuclear spin- 
lattice relaxation rate, of plane ^^Cu has the 
main contribution from the AF wavevector (tt, tt), and 
the of plane ^^O, has the contributions from 
the wavevectors in the vicinity of (tt, tt) and small g ~ 0. 
It is shown that the theory is able to explain the main 
features of experimental data on temperature and doping 
dependence of copper nuclear spin-lattice relaxation rate 
in both carrier free La2Cu04 and doped La2-a;Sra;Cu04 
compounds. 
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APPENDIX A 

Various types of sums over lattice sites utilized for cal- 
culation of frequency moments (w^) and {oj^) in the 
t—J model. The Gij and Dij symbols refer to hop- 
ping t and(or) exchange J between the first neighbors, 
whereas M and Q refer to spin-spin correlation functions 
c and(or) transfer amplitudes T at the appropriate posi- 
tions, dij is the delta function. The sums for M and Q 
were calculated for up to n = \j — i|-th neighbors (n = 
2). 

7k = {cos k^a + cos kya)/2. 

ijlmr 



V. CONCLUSION 

The theory for relaxation function in the two- 
dimensional t — J model in the paramagnetic state is 
presented taking into account the hole subsystem as well 
as both the electron and antifcrromagnetic correlations. 
The presentation obeys rotational symmetry of the spin 
correlation functions and is valid for all wave vectors 
through the Brillouin zone. The fit of effective correlation 
length S^eff to experimental data is in agreement with the 
microsegregation hypothesis, where the effect of doped 
holes affects the value of correlation length at T ^ 
as, ^0 = 1/na;, where n is the average distance between 
the holes along the domain walls. The best fit yields 



1 ^ 5„,<5i,GyD2.^yg.k.(R„-R.) ^ leMiD^G, 
ijlmr 

ijlmr 

Y SmrGuGlrGrjMlrMije''^^^"^-'^^^ = 36MfG\ 
ijlmr 

^ Y SmjGuGlmGmrMimMirf^'^^^-'-^^^ = SGMfC^lk, 
ijlmr 



ijlmr 
ijlmr 
ijlmr 

1 J2 ^ir^„^rA'rGr,•^^,•e*•(^'"-«-'■) = 16M,D^G, 

ijlmr 

^ ^ SmjSilDf^GmrMmre'''<^--^^^ = 16M,D^G^^, 
ijlmr 

^ ^ SirGlmGmrGrjMrjMime'''^^--'''-^ = MMfG^^^, 

ijlmr 

1 SirGirGimGriMi^Mrjf^'^^^--^^^ = GAMfG^'^l 

ijlmr 

ijlmr 

^ 5^ SljGimGmjGrjMimMrje'''^'^"^-'^^^ = GAMfG^^l 

ijlmr 

^ J2 SljSirGmrD^,jQmje'''<^"^-^'-^ =4D^G{3Q2 + Qo)lk, 

ijlmr 

^ 53 6mrS^JGl,Dl^Ql,e''^<^--'^^^ =AD^G{3Q2 + Qo), 

ijlmr 

1 ^ 5„,G,,Gi,G,,Mi,M,,e''^-(^--»^) = 16M,G^ (3M2 + Mq) , 

ijlmr 

^ SirGimGmrGrjMmjMire'''<^'--^^^ = 4G^ {3M2 + Mof JU, 

ijlmr 

^ Y '5yGi„G„,a,M„,Mj,e''^-(^'"-»'-) = (3M2 + Mq) (M2 (47^ - l) + Mq) , 

ijlmr 

+ 4 (QoMo + 3Q0M2 + 3M0Q2 - 7M2O2) G^7k, 



ijlmr 

^ Yl SljGijGmjGrjMimMrje'''<^'--^^^ = 16M, {3M2 + Mo)G^jI, 
ijlmr 

^ SrnrGuGirGrjMirMij^'^-^'^^-^^'^ = 4 (3M2 + Mof G\ 

ijlmr 

^ SmjGuGlmGmrMimMlre'''<^-'-^'-^ =A{3M2 + MofG^Jk, 

ijlmr 

^ Y SlrGirGmrGrjMirMmje'''^^"'-'^'''^ = I6M1 (SMa + Mq) G^k, 

ijlmr 

^ Y SijGijGimGrjM.jM^re'''-^^'--^^^ = 36M2G37k, 

ijlmr 

^ ^ SmjGi^Gi^G^rMirMi^f^'^-^^--^^ = IQM^G^ {'iM2 + Mo)^, 

ijlmr 

Jj: ^ SirGirGimGrjMirMmje'^'^^'^~^^^ = (l2 + 167^ + 8 COS fc^a COS /^ytt) , 

^ Y SijGimG„,jGrjMirMmje'''<^^-'^-'^ = M^G^ (12 + 167^ + 8 cos /c^ a cos A;^ a) , 

ijlmr 

^ Y SijGijGmjGrjMijMmre''^^^'^-^^^ = 16M,G^ (Af2 (47^ - 1) + Mo) , 

ijlmr 

1 Y SijGijGimGrjMmjMire''^^^--^'-^ = 64MlG^ji+A (M^ + 6M0M2 - 7M|) G37k, 

1 Y SmjGimGimGmrMuM^re'^<^'--'^'-^ = I6M1 (3M2 + Mq) G37k, 
1 SirGlrGmrGrjMimMrje^^<^^-'^-^ = I6M1 (3M2 + Mq) G37k, 

1 Y SlrdijDmjGrjDrjQmre'^<^'--^'-^ =AD^G{Q2{H-l)+Qo), 

ijlmr 



APPENDIX B 

In this appendix the relations to compare the result for 
frequency moments of -F(k, ui) in carrier free 2DHAF are 
given. 

For the second and fourth moments the result was^^ 

K)2DHAF = -^^f^E7.5(q), (85) 



and 



128J-^ 



K)2DHAF = -^^j:SiP)S{ci)F{Kp,cO, (86) 

where 



7p cosp^a cospj^a = 0, ^ 7q7k-q = 0, 

q 



p 



E 



7q7k-q = ^ f ^ + ^ cos /ca^a cos kya + 7^ 



E '^q'^k-q = ^"^k , E "^q 1v"'=\ 
q q 



7k-q COS qxa cos qya = ^ cos k^a cos kyU, 



F{k, p, q) ~ 7k (-27p7q - 107p7q - 47p7q + 7p7q) 
+ 7k (47p7q + 37p7q) - 27k7p7q + 57p7q 

- 7p7q + 7p7q + 57p7q7k-q + 7k+p7q 
+ 27p7k-q - 27p7k-q - 7k+p7k-q7k+p-q. (87) 

Using the relation between correlation functions in mo- 
mentum and site representation: 

S{k) = {S^Sli,) = {S^S])e'^<^'-'^^'^ = 

\i-j\<2 

= Co + 4ci7k + 4c2(47k — cos k^a cos kya — 1), (88) 

performing the summation over k and q in (86) and using 
the necessary types of integrals (sums) over the Brillouin 
zone as given below, we obtain the same expression, as 
shown in (68), with t^ff = and settling to imity the 
decoupling parameter C = 1- The integrals (sums) were 
calculated analytically and checked numerically. 



7q = (cosgra;a + cosqya)/2. 
E"''p^g4' E'>'p~I' E 



7k-q 



E 7k+p7k-q7k+p-q = Y^7k, 



2 9 

2^ 7p7k+p7k-q7k+p-q = ^7k, 

p,q 



4' 



^C0SPa;aC0SPya7k+p7k-q7k+p-q = -^7k 



2^ 7p7q7k+p7k-q7k+p-q = ^ (^k + ^7k COS fc^^o COS /cj^fl + ^7k 



45^ 



^ cospa,a cospj,a7q7k+p7k-q7k+p-q = ^ (2 COS k^a COS kya + 1) , 
p,q 



^ cospa;a cospya cos q^a cos g3^a7k+p7k-q7k+p-q = cos k^a cos A;3;a7k, 



E^p = 0' E^i^k-q = 0, ^7^=0, 
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